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ABSTRACT 

Aims. We use recent very extended HI kinematics (out to 48 kpc) along with previous Hor kinematics of the spiral galaxy NGC 3198 
in order to derive its distrihution of dark matter (DM). 

Methods. First, we used a chi-square method to model the rotation curve (RC) of this galaxy in terms of different profiles of its 
DM distribution: the universal rotation curve (URC) mass model (stellar disk + Burkert halo + gaseous disk), the NFW mass model 
(stellar disk + NFW halo + gaseous disk) and the baryon ACDM mass model (stellar disk + NFW halo modified by baryonic physics 
-I- gaseous disk). Second, to derive the DM halo density distribution, we applied a new method that does not require a global and often 
uncertain mass modelling. 

Results. While according to the standard method, both URC and NFW mass models can account for the RC, the new method instead 
leads to a density profile that is sharply disagrees with the dark halo density distribution predicted within the Lambda cold dark matter 
(ACDM) scenario. We find that the effects of baryonic physics modify the original ACDM halo densities in such a way that the 
resulting profile is more compatible with the DM density of NGC 3198 derived using our new method. However, at large distances, r 
~ 25 kpc, also this modified baryon ACDM halo profile appears to create a tension with the derived DM halo density. 

Key words. Dark Matter; Galaxy: NGC 3198; NFW haloes; Universal Rotation Curve; Baryonic Feedback 


1. Introduction 

It has been known for several decades that the kinematics of disk 
galaxies leads to a mass discrepancy (e.g. Bosma| 1 978[ |Bosma| 
& van der Kruit|[l979[ [Rubin et al.||1980| ). While in their in¬ 
ner regions that range between one and three disk exponential 
scale lengths according to the galaxy luminosity ( jSalucci & Per-j 
sic|1999| ), the observed baryonic matter accounts for the rotation 
curves (RCs) (e.g. Athanassoula et aL]|1987[ [Persic & Saluccij 
1988| [Palunas & Williams|[2000[ l, we must add an extra mass 
component in the outer regions, namely a dark matter (DM) halo 
to account for that component. The kinematics of spirals is now 
routinely interpreted in the framework of a DM component. In 
the widely accepted Lambda cold dark matter (ACDM) scenario, 
the virialized structures are distributed according the well known 
NFW profile proposed by Navarro, Frenk, and White ( [Navarro [ 
et al. 1996| l. The ACDM scenario describes the large-scale struc¬ 
ture of the Universe well (e.g. [Springel et al.|2006|l, but it seems 
to fail on the scales of galaxies ( de Blok & Bosm^[2002[[Gen-[ 
tile et al.|2004 2005 [l. Going into detail, the NFW density pro¬ 


file leads to the “core-cusp problem”: empirical profiles with a 
central core of constant density, such as the pseudo-isothermal 
([Begeman et al. 1991[ Kent[[1986 l and the Burkert (Salucci & 


[Burkert[[2000[ ), fit the available RCs much better than the mass 
models based on NFW haloes. 

In the present paper, we derive the DM content and distri¬ 
bution in the spiral galaxy NGC 3198. This galaxy has been the 
subject of several investigations. It was studied by means of opti- 


and HI-21 cm radio observations ( |Bosma| 1981 [[van Albada et al.[ 
[1985[[Begeman|1987] the latter established it as the object with 
the clearest evidence for DM), see also ([de Blok et al. [[20081 
[Gentile|2008| l. 

Our present analysis is mainly based on the HI observations 
by [Gentile et al.[p013[ ), part of the HALOGAS (Westerbork Hy¬ 
drogen Accretion in LOcal GAlaxieS) survey. The main goal 
of HALOGAS is to investigate the amount and properties of 
extra-planar gas by using very deep HI observations. In fact, for 
this galaxy, they present a very extended RC out to 720 arc- 
sec, corresponding to ~ 48 kpc for a galaxy distance at 13.8 
Mpc ([Freedman et al.[[200r| l. The previous HI observations by 
de Blok et al. (20081 were only extended out to ~ 38 kpc, for 
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the same galaxy distance. In [Gentile et al.[ ( (2013[ ) this extended 
RC was modelled in the framework of modified Newtonian dy¬ 
namics (MOND). Here, we want to use such a uniquely extended 
kinematics to help resolving the DM core-cusp issue. To the very 
reliable kinematics available from 2 to 48 kpc we apply two dif¬ 
ferent mass decomposition methods that will derive the DM halo 
structure. This is compared with a) the empirically based halo 
profiles coming from the URC, b) the NFW haloes and c) the 
baryon ACDM haloes, the outcome of scenarios in which bary¬ 
onic physics has shaped the DM halo density. 

This paper is organized as follows. In Sect. 2 we present the 
HI and Ha kinematics used in this study. In Sect. 3 we model 
the RC by using the quadrature sum of the contributions of the 
individual mass components (stellar disk -H dark halo -H gas disk) 
where the dark halo has a NFW or a Burkert density profile, 
respectively. In Sect. 4 we obtain the results of standard mass 
modelling of the NGC 3198 RC. In Sect. 5 we use a recent mass 
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Fig. 1. Comparison between Hu and HI RCs black open triangles with error bars from Corradi et al. (1991 
[Daigle et al.|(|2006t. and red circles with error bars are fromjGentile et al.](|2013^. 


blue circles with error bars are from 


modelling technique to obtain a very robust and careful determi¬ 
nation of the DM halo density of NGC 3198 and to show that it 
is at variance with the NFW density prohle in an unprecedented 
way. We consider the mass dependent density prohles obtained 
by hydro-dynamical simulations of ACDM haloes in Sect. 6. Our 
conclusions are drawn in Sect. 7. 


2. KINEMATICS DATA 

2.1. HI data 


The HI data of NGC 3198 were taken in the framework of the 
HALOGAS sur vey ([Heald et al.|2011 1 and they were presented 
in jGentile et al.| ( |2013| l. The data were obtained with the WSRT 
(Westerbork Synthesis Radio Telescope) for 10x12 hours, with 
a total bandwidth of 10 MHz subdivided into 1024 channels. The 
data cube used to derive the rotation curve has a beam size of 
35.2 X 33.5 arcsec, and we were able to detect emission down to 
~ 10'^ atoms cm“^. 

To construct the gas distribution and the rotation curve of 
this galaxy in a reliable way , we modelled the whole data cube 
by means of the TiRiFiC software (Jozsa et al. 2007). We suc¬ 
cessfully modelled the HI observations of NGC 3198 with a thin 
and a thick neutral hydrogen disk, and, thanks to an increased 
sensitivity, we were able to trace the rotation curve out to a dis¬ 
tance of ~ 48 kpc (for the galaxy distance of 13.8 Mpc), i.e. to a 
larger radius than those reached in previous studies. More details 
about the data reduction analysis and modelling can be found in 
[Gentile et al.| ( |2OT3l l. 


2.2. Ha data 


The Ha rotation curves of NGC 3198 have been published by 


many authors (jCheriguene |1975| Hunter et aL||1986( 

Bottema 

19881 Wevers et al.|1986||Kent|1987| Corradi et al.|199 

[[[Daigle 

et al.|2006 1. We notice that Corradi et al.|([1991[l and[Daigle et al. 


( |2006|l measurements are a good representation of these data. In 
Fig. [ij these RCs are plotted along with the HI RC used in this 
work. Ha data provide us measurements of the circular veloc- 
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Table 1. Stellar disk contribution Vj (km s“') ([de^loketalJ2008jl and 
the circular velocity V (km s“*) mainly from Gentile et al.| ( j201^ , but 
also from [Corradi et ar] l [1991[ > and [Daigle et al.| ( [2006^ , of the NGC 
3198 with errors dV (km s ')' The DM density prohle p x 10“^^ (g 
cm“^) (see Sect. 5 for details). 


R 

(kpc) 

Vd 

(km s-') 

V 

(km s-') 

dV 

(km s-') 

P X 10-25 

(g cm-5) 

2.0 

86.2 

79.0 

7.0 


3.0 

85.4 

97.8 

5.0 

— 

4.0 

93.6 

118.0 

5.6 

— 

5.5 

115.7 

139.4 

4.3 

2.34 

6.0 

120.8 

144.2 

4.3 

2.33 

7.0 

125.4 

143.3 

4.5 

2.20 

8.0 

125.5 

150.3 

4.3 

2.01 

9.0 

123.5 

149.9 

4.3 

1.83 

10.1 

120.1 

152.1 

4.3 

1.64 

11.0 

116.6 

151.1 

4.5 

1.48 

12.1 

112.6 

156.2 

4.3 

1.32 

14.1 

105.2 

161.0 

4.3 

1.06 

16.1 

98.6 

155.3 

4.3 

0.86 

18.1 

92.7 

148.7 

4.3 

0.70 

20.1 

87.5 

149.1 

4.3 

0.58 

22.1 

82.8 

148.4 

4.3 

0.48 

24.1 

78.7 

146.2 

4.3 

0.42 

26.1 

75.1 

145.5 

4.3 

0.36 

28.1 

71.9 

147.3 

4.3 

0.33 

30.2 

68.9 

146.5 

4.3 

0.30 

32.2 

66.3 

148.4 

4.3 

0.27 

34.2 

63.9 

149.3 

5.0 

0.25 

36.2 

61.8 

149.9 

4.3 

0.23 

38.2 

59.8 

149.3 

4.3 

0.21 

40.2 

58.0 

150.0 

4.6 

0.20 

42.1 

56.4 

147.6 

7.0 

0.18 

44.2 

54.9 

149.8 

4.3 

0.16 

46.2 

53.5 

151.5 

4.3 

0.13 

48.2 

52.2 

151.9 

7.7 

0.11 


10 
















































E.V. Karukes et al.; The dark matter distribution in the spiral NGC 3198 out to 0.22 R^ir 


ity at five different new (inner) radii not mapped by the present 
HI RC and four more measurements at radii at which we can 
combine them with our HI data. The hybrid RC (HI+Ha data) is 
listed in Table [T] All together, HI and Ha RCs agree very well 
within the observational errors where they coexist. With respect 
to the HI RC of |Gentile ( |2008 1 , the present hybrid RC covers the 
innermost and the outermost regions significantly better; how¬ 
ever, Ha RCs give us no new useful information for r < 2 kpc, 
owing to its large observational uncertainties and because in this 
very inner region the kinematics is strongly influenced by non- 
axisymmetric motions (see Corradi et al.| 199 i'll. 


3. MASS MODELS 


3.2.1. Burkert halo 


The URC of galaxies and the kinematics of individual spirals 
(|Salucci et al.|2007)l points to dark haloes with a central constant 


density core, in particular, to the Burkert halo profile (Burkert 
1995[|Salucci & Burkert|2000) l. The relative density distribution 


is given by 


PuRc(0 - 


Pore 

(r + rc)(r2 h- r^) 


( 2 ) 


where po (the central density) and (the core radius) are the two 
free parameters. The present data cannot distinguish these URC 
haloes from other cored profiles, for which lim^^o p{r) = const. 


We model the spiral galaxy NGC 3198 as consisting of three 
“luminous” components, namely the bulge and the stellar and 
the gaseous disks, which are embedded in a dark halo. To study 
the properties of luminous and dark matter in this galaxy, we 
model the RC in terms of the contributions from the stellar disk, 
the bulge, the gaseous disk, and the dark matter halo: 


3.2.2. Navarro, Frank and White (NFW) halo 

In numerical simulations performed in the (A) CDM scenario of 
structure formation, [Navarro et al.| ( |1996] l found that virialized 
systems follow a universal DM halo profile. This is written as 


vHr) = Vlir) + Vlir) + V^ir) + y2^(r). 


( 1 ) 


3.1. Luminous mailer 


We define Vd(r) as the contribution of the stellar disk to the cir¬ 
cular velocity V(r). The surface brightness profile of NGC 3198 
has been analysed very well by jde Blok et al.| ( [2008) l. We assume 
their one-component surface brightness profile to derive Vd(r). 
As a reference value w e take the stellar mass-to-light ratio from 


de Blok et al. 




= 0.8 (referred to the Spitzer 
IRAC i.b pm band, which is a good proxy for the emission of 


the stellar disc). We thus set V^ir) - (Vj®®’°’^(r))^-QY- Then, we 
leave the amplitude of the disk contribution to the circular veloc¬ 
ity (i.e. the disk mass) as a free parameter to be derived by fitting 
the RC. No results of the present paper will change by assuming 
any other models of the dostribution of the stellar disk of NGC 
3198 in jde Blok et al.| ( |2008| l or in previous works. 

The contribution from the bulge component is Vh(r). We fol¬ 
low jde Blok et al.j ( |2008| ) and we consider their 1 -component 
model which accounts for this inner stellar component. Hence, 
Vd{r) - 0 since the bulge is included in y^®®*°^(r). 

The helium corrected contribution of the gaseous disk de¬ 
rived from the HI surface density distribution V^{r) taken from 
[Gentile et al.| ( |2013| l. We notice that, thanks to the accuracy of 
the HI measurements and to the excellent knowledge of the dis¬ 
tance of this galaxy, this component is derived here much better 
than in the majority of the spirals studied to resolve the core-cusp 
controversy. 

The "luminous" component of the circular velocity of this 
galaxy is then well known, except for the value of the disk mass 
which contributes to putting NGC 3198 at the front line of DM 
research. 


PNFw( 0 - _r 

where ps and r^ are the characteristic density and the scale ra¬ 
dius of the distribution, respectively. These two parameters can 
be expressed in terms of the virial mass Myii, the concentra¬ 
tion parameter c - and the critical density of the Universe 
Peril = 9.3 X 10“^°gcm“^. By using Eq. we can write 


P.9 

rs 


100 

3 log(l-Hc)- 


1 / 3 X , 

C \47rl00pcritj 


gem 


(4) 


where c and are not independent. It is well known from 
simulations that a c - relationship emerges ([Klypin et al.[ 
[20TT] [Bullock et al.|2001[|Wechsler et al.|2002| l: 


11.7 




10"Mp 


- 0.075 


(5) 


We are testing the density profile of haloes made by 
collision-less cold dark matter particles. Variations in this sce¬ 
nario are not considered in this work, except one in Sect. 6. 


4. RESULTS FROM THEFITTING METHOD 

The mass modelling results for the Burkert and NFW profiles are 
shown in Figs (|^ and ([^, respectively. The URC Burkert halo 
gives an excellent fit (see Fig.|^ with a reduced chi-square value 
ofx^ - 0.9. The best-fit parameters are 


3.2. Dark mailer 

r 

We define yDM(0 = f 47rpDMR^dR as the contribution to 
0 

y(r) from the dark matter halo of the virial mass M,,,> = 
jTrlOOpcri/Ryjy that could have a variety of density profiles. Here 
Pdm represents a DM halo density profile. 


po = (3.19 + 0.62) X 10^^^ gem^^ 
re = (17.7 + 2.0) kpc; 

= (0.98 + 0.04). 


Then, we compute the mass of the stellar disk as 


Md - l.lT^®(yd)(Ri)2| 


(6) 
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Fig. 2. URC mass modelling of NGC 3198. Circular velocity data (filled circles with error bars) are modelled (thick red line) by the halo cored 
component (thick green line), the stellar disk (magenta line) and the HI disk (azure line). 



Fig. 3. NFW mass modelling of NGC 3198. Circular velocity data (filled circles with error bars) are modelled (thick red line) by the stellar disk 
(magenta line), the NFW halo profile (green line) and the HI contribution (azure line). 


where VdC^i) is the disk contribution to the circular velocity at 
the outermost radius R] 48 kpc. This estimate is very solid 
and independent of the actual light profile in the inner part of the 
galaxy. We find Md - 4.4 x IO'^Mq (=s 2 times bigger than the 
values found by |de Blok et al.||20()8] l with a propagated uncer¬ 
tainty of about ten percent. The corresponding virial mass and 
virial radius are Mvh- = 5.8;|;Qg x 10"Mo and 
kpc.The 1,2,3-cr confidence regions for the best-fit parameters 
are shown in Fig. The central points correspond to the best-fit 
values. 

In the framework of the NFW mass models, we fitted data 
in terms of the free parameters; the virial mass, the concen¬ 
tration parameter and above defined the mass-to-light ratio 
(Mvir, c, TJ ®). The results of the best-fit are 
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= (8.9 + 2.1) X 10" Mo; 
c = (6.69+ 1.46); 

= (0.79 ± 0.07). 

In this case the reduced chi-square is - 0.8, even slightly 
better value than found for the URC-halo model. The best-fit 
value of the concentration parameter c - 6.69 ± 1.46 is found 
to be somewhat lower than what is expected from Eq. 
cnfw ~ 10 + 1. It is worth recalling that in other galaxies, this 
discrepancy in the concentration parameter is much larger (see 
|McGaugh et al.]2003[|Salucci et al.|2010||Memola et al.|2011) l. 

From Eq. the disk mass, within a ten percent uncertainty, 
is Md = 3.5 x IO'^Mq, a somewhat smaller value than found 
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Fig. 4. 1,2,3 cr confidence ellipses (purple, red, orange, respectively) of the best-fit parameters in the Burkert halo case. The central points indicate 
the best-fitting values. “ is in the IRAC 3.6 pm band units. 



Fig. 5. 1,2, 3cr confidence ellipses (purple, red, orange, respectively) of the best-fit parameters in the NFW halo case. The central points indicate 
the best-fitting values. “ is in the IRAC 3.6 pm band units. 


for the URC-halo model. The best-fit values of the scale radius 
and the characteristic density are Tj = (37.2 + 11.0) kpc and 
Ps = (8.0 + 4.1) X 10“^® g cm~^. 

The 1,2,3-cr best-fit parameters confidence regions are 
shown in Fig. The central points correspond to the best-fit 
values that result somewhat higher or lower than the N-Body 
simulation outcome relative to a galaxy with V^ax - 150 km/s 
as NGC 3198. The discrepancy, however, is within 1.5 cr. 

The standard mass modelling of the kinematical data of NGC 
3198 is then not able to clearly select between a cored and a 
cuspy halo profile. In fact, in the case where a galaxy like NGC 
3198 showing a flattish rotation curve over a wide range of radii, 
we have a modelling degeneracy; the same best-fit solution cor¬ 
responds to very different mass models (see Appendix of|Gentile| 
let al.|20041 l. 


where a^i, and ahi the radial accelerations, generated, re¬ 
spectively, by the halo, stellar disk, and HI disk mass distribu¬ 
tions. Within the approximation of spherical DM halo, we have 


flH = 47rGr ^ J"pH{r)r^ dr. 
0 


( 8 ) 


Therefore, we compute the derivative of Eq. (0 by manipulating 
previous equations to get the DM density at any radius. We have 


Pnir) = 


Xq d 
AnGr^ dr 




(9) 


5. A NEW METHOD OF ESTIMATING THE HALO DM 
DENSITY AND ITS RESULTS 

A step forward in mass modelling spirals has come from the 
method of |Salucci et al.|f2010) l. This method was applied first 
to the Milky Way to derive the value of the DM density at the 
Sun’s location. In this paper we are applying it, for the first time, 
to the outermost parts of an external galaxy. Our aim is to derive, 
for the spiral galaxy with the most extended kinematics, the DM 
density at large radii, where the influence of the stellar and HI 
disks is respectively negligible and known. 

The idea underlying the |Salucci et al.| ( |2010| l method is to 
resort to the equation of centrifugal equilibrium holding in spiral 
galaxies: 


y2 

- — Ah + 

r 


(7) 


where is a factor correcting the spherical Gauss law used 
above in case of any oblateness of the DM halo. Since this value 
is very near to one, we assume - 1 (see details in jSalucci 
let al.|20101 i. 

Equation(j^ gives a very good estimation of the density when 
the contribution from the luminous components is small. In 
short, to work well the new method requires a high-resolution, 
high-quality, and very extended kinematics for a spiral of known 
distance. Of course, we are interested in the region well outside 
3Rd. 

Eor simplicity, we model the disk component as a Eree- 
man stellar exponential, infinitesimally thin disk ( |Ereeman|1970| l 
with the disk scale length Rr, - 3.7 kpc. No result in this paper 
will change by instead assuming the de Blok et al. ( |2008| l disk 
mass profile, whose contribution to the circular velocity Vd is 
given in TablefT] Eor Mq we use the URC mass modelling value: 
Md - 4.4x IO^Mq, and no result changes if assuming any other 
reasonable value for this quantity (see below). 
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Fig. 6. Logarithmic slope of the RC of NGC 3198. 


The surface stellar density profile is given by 


2(r) = 


^ Mp ^ 
ylirRlj 


e "d . 


We can write 


( 10 ) 


GMor 

aD(f)= —^(/o^o-/i^fi) 
"d 


( 11 ) 


where /„ and are the modified Bessel functions computed at 
rl(2R^). 

The HI disk componen t Vg(r) and its deri vative are easily 
obtained from observations (Gentile et al. 20131. Finally, Eq. (j^ 
becomes 






Ri 


\Rd 


—^(1 + 2y(r))J 


( 12 ) 


where = (3/o.^fo - ^i^i) + ^(A^o -^o^i) and a{r) and 

y{r) are the logarithmic slopes of the circular velocity and of the 
HI+He disk contribution to the latter, both of which known. 


We stress that in galaxies with VdORi)) ^ y(3/?D)j it 
is very difficult to use the standard mass modelling method 
to disentangle the circular velocity into its dark and luminous 
components and to obtain the DM density distribution out to 
r ^ 6/?d- Instead, the fundamental point of the present new 
method is that, for radii r > 3/?d the second term of RHS of 
Eq. ( [T2) l always goes rapidly to zero becoming much smaller 
than the first and the third terms, both known. Then, by means 
of Eq. ( [T2j l, we can immediately derive p{r): the unknown 
term, proportional to the stellar disk mass becomes irrelevant 
as r > 3 Rd- In Eq. ( [T2| , all terms have the dimensions of a 
density; specifically, the three terms of RHS can be considered 
as the effective density of the whole gravitating matter and the 
(sphericized) densities of the stellar and gaseous disks. 

Thus, by means of Eq. 0 for r > 3 Rd, we obtain a re¬ 
liable dark matter density profile, the farther we get, the more 
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precise the estimate becomes. In the range : r < 3Ru, Eq. 0 
still holds, but it cannot give information about the DM halo 
since this is very sub-dominant there. In this region, however, 
we can use the latter equation to derive the stellar disk mass. 
This estimate, however, may turn out to be somewhat uncertain 
because, in this inner region, the first term of the RHS of Eq. 
( [T2] ) has some observational errors; moreover, moderate errors 
in the measured value of the disk scale length can affect the third 
term. Incidentally, we notice that the distance of the galaxy must 
be known with good precision (as it is in NGC3198), because its 
uncertainty affects all three terms of the RHS in different ways. 

In the case of NGC 3198 we take y(r) = 0 for r > 3 Rd, in 
order to simplify our calculations, since the gas contribution to 
the circular velocity, from ~ 12 kpc, is nearly constant (see Eig. 
1^. The logarithmic slope of the RC instead varies with radius 
(see Table[2l, even at r ~ 17 kpc: we take dlog V(r)/dlogr as in 
Eig. 1^ In Eig. |7] we obtain the density profile of the dark halo. 
We list the values of the obtained DM profile starting from 5.5 
kpc in Table [2 We see that the DM component starts to dominate 
the luminous components from r ~ 10 kpc; moreover, the stellar 
disk’s contribution in the RHS of Eq. goes further below 
the gravitating matter for r > 10 kpc and to zero for r > 17 
kpc, independently of its mass. This means that, starting from 
~ 17 kpc, the halo density profile, which obtained by means 
of Eq. ( [T2 ] i, is virtually free from the uncertainty on the actual 
value ofthe disk mass, which usually plagues the standard mass 
modelling of RCs. 

Next we discuss in detail the distribution of matter in the var¬ 
ious regions of NGC 3198. In the innermost one, r <2 kpc we do 
not have any kinematic al information due to the lack of data. In 
the region extended from ~ 2 kpc to ~ 10 kpc, the stellar compo¬ 
nent dominates over the DM component. No direct information 
on the latter can be extracted here. In this region, however, we 
can use the fair agreement between the stellar and the dynamical 
density (as defined in Eq. ( [T^ ) to derive the value of Mo- 

Eor r > 17 kpc, the DM density is directly obtained by means 
of Eq. ( |T2 ) i see Eig. Here, all quantities entering in Eq. ( [T2) i 
are known within a small uncertainty, and this leads to a robust 
determination of p{r). 

The effect of the uncertainty in the value of the disk mass 
in the determination of p{r) is shown well in Fig. We now do 
not take a specific value for Mo, and we derive the DM density 
profiles through Eq. ( [T^ by assuming different values for the 
disk mass. We see that the |Salucci et al.| ( f2010[ ) method leads to 
a family of density profiles that all agree outside the stellar disk 
(i.e. for r > 3Rd), independently of the corresponding value of 
the disk mass (see the zoomed area in Fig. . 

We compare the derived DM density with the URC and the 
NFW profiles, see Fig. In the first case, the derived density 
bears no difference with that obtained by means of the first 
method, i.e. with Eqs.([T]) and 0. We therefore found Mo = 
(4.4 + 1.0) X IO^Mq. This is also evident in the zoomed area 
of Fig. [8| (see Fig. [9] ). In external regions of the NGC 3198, the 
[Salucci et al.| ( |2010) l method yields a halo that is compatible with 
the URC halo, as derived by chi-square fitting of the RC of NGC 
3198. 

For NFW haloes the situation is very different. First, in Fig. 
I^we realize that, independently of the disk mass we assume, the 
best-fit NFW halo profile is in poor agreement with the derived 
density. We consider this fault only as a hint and not as a defini¬ 
tive evidence against a NFW halo. Within the new modelling 
method, we derive the value of c and Myir and their uncertainties 
by evaluating, from Eqs. 0, (0, and ( [T^ , the DM density (and 
its uncertainties) at the radius r = 45 kpc (where the stellar disk 
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Fig. 7. Density profile of the DM halo of NGC 3198 and the effective density of the other components. We assume Mq = 4.4 x IO’^Mq. The 
stellar disk (Blue line), the HI disk (magenta line), the dark halo (red line) and the sum of all components (green line). The different colour regions 
correspond to the regions a) where we do not have any kinematical information due to the lack of data (dark purple), b) where the stellar disk 
dominates the DM density profile (light purple) and c) where DM dominates (white). 



r(kpc) 

Fig. 8. DM density profiles for different disk mass values Md = nx IO'^Mq and n=2,3,4.4,5 from the highest to the lowest curve (green, blue, azure 
and magenta triangles). The black solid line shows the NFW density profile with Myij = 8.9 x 10 '*Mq corresponding to the best fit values found 
in the previous section. Two black dashed lines show the NFW density profiles taking Icr uncertainties in Mv,t and 10 percent of uncertainties in c 
into account. The URC halo is shown as Red solid line. Two Red dashed lines show the URC density profiles obtained by taking Icr uncertainties 
in po, Cc into account. 


certainly does not contribute to the gravitating matter density 
profile). A more serious problem appears when we realize that 
the resulting values of the halo parameters and c are very 
different from the ones derived by means of the standard method 
applied in the previous section. In fact, for the virial mass and the 
concentration parameters we get Myir = (8.9 + 2.1) x 10’'Mo, 
c = 6.69+ 1.46. 

In Fig. [T^we plot the two solutions for (c - Myij) with their 
uncertainties. We now compare these values of the concentration 
parameter and the virial mass (c - Myu) obtained by each of the 


mass modelling techniques also in light of the {c-Myff) relation¬ 
ship that emerged from numerical simulation (see Eq. (|^). They 
agree only for values of c that are much lower than those emerg¬ 
ing in numerical simulations and for values of My^ > 9xlO"Mo 
far too high for this spiral galaxy which has y^ax ^ 150 km/s. 

In short, by assuming a NFW halo in NGC 3198, we have 
two completely different solutions for its structure parameters, 
according to whether one adopts the local or the global method 
of mass modelling. The best intersection of these sets of re¬ 
sults is in total disagreement with the properties of the simulated 
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Fig. 9. Zoom of Fig.[^on linear scales. 



Fig. 10. NFW case: the 1,2,3 - cr (purple region, red region and orange region, respectively) confidence ellipses for the global best-fit parameters. 
The blue solid line shows the c - relation from numerical simulations, the dashed blue lines show its 10 percent uncertainty. The Green region 
shows the c — My^ relation from the local density values obtained by the|Salucci et al.|j2010^ method taking the 10 percent uncertainty into account. 


haloes in the N-body ACDM scenario. Noticeably, the problem 
for NFW haloes is different here and, if possible, more seri¬ 
ous than that of the core-cusp discrepancy, usually occurring at 
~ 0.05f?vir ( |Donato et al.|200^ . We found, in fact, that the den¬ 
sity of the DM halo around NGC 3198 is not very consistent with 
the NFW profiles well out to 0.22/?vir. 

Of course, NFW haloes emerge out of DM only simulations 
in the leading ACDM scenario. However, this scenario actual 
haloes around galaxies may have undergone a compression by 
the stellar disk during the formation of the latter (e.g. Gnedin | 
et al. 2004[ Katz et al. 2014|l and/or suffered a baryonic feed¬ 


back during the subsequent history of the galaxy ( |Di Cintio et al.| 
|2014| l. These processes could have modified the DM haloes orig¬ 
inal distribution. Furthermore, haloes today around galaxies, can 
be born within a different cosmological scenario like Warm or 
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Self-Interacting DM, in which a NFW profile is not established 
(e.g. |de Vega et al.|2014[ ). 

In spite of that, in this section we focused on the failure of the 
NFW halo profile because: a) it is important to fully expose the 
discrepancy between the density distribution of the dark haloes 
around galaxies and the predictions of the simplest dark particle 
scenario; b) NFW haloes are still used often to investigate impor¬ 
tant Cosmological issues in the belief that the discrepancy with 
actual galaxy haloes, though present, does not have much phys¬ 
ical relevance; c) de facto, several cosmological investigations 
have been carried out considering that there is no discrepancy at 
all. 

It goes without saying that the failure of NFW profiles is not 
the demise of the ACDM scenario (see next section). 
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Fig. 11. Observed DM halo density profile (black points) of NGC 3198 
ICintio et al.H2Q14| mass-dependent model predictions (red line). 


6. Testing LCDM halo profiles modified by the 
physics of stellar disk formation 

A line of thought holds that the cosmological core-cusp problem 
(CCCP) can be addressed by considering the present DM haloes 
around galaxies like NGC 3198 as very different with respect to 
those emerging out of N-body simulations. In fact, the forma¬ 
tion and the growth in them of stellar disks and the related nu¬ 
merous and powerful supernova explosions could have modified 
the original N-body profile, by making it shallower and more in 
agreement with observations. We stress that such a dissolution 
of the CCCP is rarely studied. 

One exception is the work by |Di Cintio et al.| ( 2014[ ) based 
on the analysis of hydro-dynamically simulated galaxies drawn 
from the MaGICC project ( [Brook et al.||2012{ [Stinson et al.j 
[2013[ l. They did find, at the end of these simulations, that DM 
haloes had a completely new family of profiles, a generic double 
power-law density profile ( [Di Cintio et al.[20T4| ). This was found 
to vary in a systematic manner in the stellar-to-halo mass ratio 
of each galaxy. 

The mass-dependent density profile was derived by 
analysing hydrodynamical cosmological simulations. This pro¬ 
file (hereinafter referred to as DC 14) accounts for the effects of 
feedback on the DM haloes due to gas outflows generated in high 
density star-forming regions during the history of the stellar disk. 
The resulting radial profile is far from simple, since it starts from 
an (a,/?, y) double power-law model (see[Di Cintio et al.|2014[l 


r(kpc) 


and the obtained errors by the error propagation analysis (purple area). Di 


a = 2.94-logio[(10^+^'^V‘®® 

P = 4.23-I-1.34Xh-0.26X^ (14) 

y = -0.06-Hlogio[(10^+2.56)-0.68^(jQX+2.56^] 


where X = logio(^). 

The concentration parameter of the halo is c = ^. An alter¬ 
native definition, adopting the radius r_ 2 , is cdci 4 = where 
r _2 is the radius at which the logarithmic density slope of the 
profile is -2. This definition allows defining the same physical 
for different values of (a,/3, y) 



(15) 


Following [Di Cintio et al.[ ( |2014[ l, we reach the relation be¬ 
tween the concentration parameter in the hydrodynamical simu¬ 
lations and the N-body haloes as a function of the stellar-to-halo 
mass ratio: 


CDC 14 = (1-0 + 0.00003e^ '*'^^'‘ ®^) x cnfw- (16) 


PdC14(^) - 


(^)"(l + (^)“)“ 


(13) 


where ps is the scale density and the scale radius. The inner 
and the outer regions have logarithmic slopes -y and -/?, re¬ 
spectively, and a indicates the sharpness of the transition. These 
three parameters are fully constrained in terms of the stellar-to- 
halo mass ratio as shown in[Di Cintio et al.[P014[l: 


Using the definition of the enclosed mass, we can write down 
the expression for the scale density of the DC 14 profile: 


Mv 


47r f 






dr 


(17) 


To reduce the number of free parameters of the DC14 profile, 
it is necessary to adopt a value for the concentration parameter 
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of the original NFW halo. For our purpose and since this quan¬ 
tity depends extremely weakly on the virial mass, we can adopt; 
cnfw = 10. By bringing together all the above equations of this 
section and noticing that = CDci 4 f- 2 , we can rewrite Eq. ( [T3] l 
just as a function of the scale radius and the stellar-to-halo mass 
ratio. 

With this density profile, we can attempt a two-free parame¬ 
ter (X, rs) fit of the derived halo density profile of the NGC 3198 
(Fig. 0- The result is in Fig.[TT] the best-fit parameters are: 


X = -2.6; 

Vs — 12.2kpc. 

The inferred virial radius is /?vir - 235kpc and virial mass 
is Mvir - 7.7 X IO^Mq. As result, the CCCP is only partially 
resolved, and the DM core we detect in NGC 3198 can be ex¬ 
plained as an effect of the energy that, in different ways, stars 
have injected in the galactic ambient. However, in the outer den¬ 
sities of the dark halo derived for the first time in this study, there 
are features that may conflict with the N-body simulation predic¬ 
tions that should be recovered at /? > 25 kpc from the centre of 
NGC 3198. 


7. Conclusions 

Galaxies with a flattish rotation curve between 5 kpc to 50 kpc 
(e.g. the spiral NGC 3198, for many years the flagship of the 
evidence of the DM in galaxies) amount to only a few percent of 
the total number of disk systems. However, they play a crucial 
role in the core-cusp issue of the DM density. In fact, these "flat" 
RCs can be well fitted by either a cored DM halo or a cuspy DM 
halo, just by adjusting the amount of the stellar matter content. 
In contrast, in spirals in which |dlog V(r)/dlogii is significantly 
far from zero, this quantity constraints the distribution of DM, 
usually towards a cored one ( [Gentile et al.|2005j l. 

NGC 3198 is a special galaxy. The HI disk of this galaxy is 
very extended out to ~ 13 disk scale lengths or out to ~ 0.22 
Rvir- Furthermore, for the innermost region, many complemen¬ 
tary optical kinematical measurements are available. This object 
(of known distance of 13.8 Mpc) shows a spectacular evidence of 
a dark force in action; baryons of this galaxy are clearly unable 
to account for its (very extended) kinematics. We must assume 
that a large part of the circular velocity of NGC 3198 is due to 
a DM halo. In addition, the circular velocity of this galaxy is at 
variance with the MOND paradigm Gentile et al. ( 2013|l, while 
it seems plausible within the F(R) scenario (see [Salucci et al.| 
|2014| l. It is then obvious that to resolve the core-cusp issue in 
NGC 3198 is of particular importance. 

To do so, we used the old optical and new HI kinematics 
performed in the HALOGAS survey in combination with a new 
method of mass modelling a rotation curve. First, we verified that 
with the standard mass modelling of the RC, both URC Burk- 
ert and the NFW dark halo models fit the available data well. 
Then, once assured of the conditions of its applicability, we used 
the new refined method developed by |Salucci et al.| ( |2010[ ) to de¬ 
termine the DM halo density of NGC 3198. This result is incon¬ 
sistent with the NFW haloes predictions, independently of any 
assumption about the luminous component we can take. More¬ 
over, the derived density profile strongly supports a cored distri¬ 
bution of DM, adding independent evidence to the idea of cored 
DM distribution in galaxies. 

Today, within the ACDM scenario, NFW haloes are often 
still assumed, although it has been recognized that in this sce¬ 
nario the actual DM halo density is not the one emerging from 
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cosmological N-body simulations. In fact, it has been agreed that 
they have been modified by the subsequent baryonic physics. 


leading to stellar disk formation and evolution. Di Cintio et al. 


P014| l have simulated this important phase of the history of spi¬ 
rals and found the emerging profiles of dark matter haloes. This 
baryonic ACDM halo profile prediction fits the detected halo of 
NGC 3198 very well, especially in its cored region. At very large 
distances, 25 kpc, however, the DM halo density derived here re¬ 
sults in a clash; i.e., it is significantly higher than the outcome of 
the hydrodynamic N-body ACDM simulations. This disagree¬ 
ment is not an isolated one (see [Gentile et al.|2007|l. 
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